function y=normeL2(t,u)
    z=u.^2;
    y=sqrt(trapz(t,z));
    %end